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Abstract  The  method  of  vector  Lyapunov  functions  to  determine  stability  in 
dynamical  systems  requires  that  the  comparison  system  be  quasimonotone 
nondecreasing  with  respect  to  a  cone  contained  in  the  nonnegative  orthant.  For  linear 
comparison  systems  in  R"  with  real  spectra,  Heikkila  solved  the  problem  for  n  =  2  and 
gave  necessary  conditions  for  n  >  2.  We  previously  showed  a  sufficient  condition  for 
n  >  2,  and  here,  for  systems  with  complex  eigenvalues,  we  give  conditions  for  which  the 
problem  reduces  to  the  nonnegative  inverse  eigenvalue  problem. 
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1.  Introduction 

The  technique  of  vector  Lyapunov  functions  has  proven  very  useful  in  analyzing  the 
stability  of  differential  systems  (see,  e.g.,  [9]).  In  this  method,  however,  the  comparison 
system  requires  stability  properties  as  well  as  quasimonotonicity  relative  to  a  cone  (see 
[8]).  For  linear  comparison  systems,  it  is  sufficient  to  consider  square  systems  A(t)  £ 
R"xn,  t  >  0,  and  to  find  a  cone  contained  in  R_p  with  respect  to  which  the  system  is 
quasimonotone  (see  [1]). 

Heikkila  [6]  showed  that  to  solve  the  problem  for  linear  comparison  systems  A  £  Rnxn, 
it  is  necessary  and  sufficient  to  find  a  nonnegative,  nonsingular  matrix  B  £  Rnx"  such 
that  if  C  =  B~1AB,  then  Cij  >  0  for  all  i  ^  j.  In  this  case  we  say  that  A  is  cross-positive 
on  the  proper,  simplicial  cone  determined  by  the  columns  of  B  (see  [10])  and  that  C  is 
essentially  nonnegative  (see  [11]). 

In  [6],  Heikkila  used  Perron-Frobenius  theory  to  show  a  necessary  condition  is  for  the 
matrix  A  to  have  a  real  eigenvalue  with  greatest  real  part  and  an  associated  nonnegative 
eigenvector.  For  n  =  2,  Heikkila  showed  that  it  is  necessary  and  sufficient  that  A  have 
two  nonnegative  eigenvectors.  For  n  >  2  he  further  showed  that  it  is  sufficient  for  a 
matrix  A  with  a  real  spectrum  to  have  all  of  its  eigenvectors  nonnegative,  and  in  [7], 
Koksal  and  Fausett  extended  this  result  to  include  generalized  eigenvectors.  In  [1],  we 
showed  that  with  n  >  2,  for  matrices  with  real  spectra  it  is  sufficient  that  this  first 
eigenvector  be  positive,  and  we  discussed  cases  when  it  is  nonnegative  in  terms  of  the 
reducibility  of  the  matrix  A. 
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Quasimon.  Of  Lin.  Sys. — Complex  Spectrum 


This  problem  has  further  applications  to  the  positive  orthant  stabilizability  problem 
(see  [3])  and  the  “hit  and  hold”  problem  (see  [2])  from  control  theory.  In  this  paper  we 
address  the  case  of  the  general  spectrum  (when  the  eigenvalues  of  A  are  not  strictly  real) 
and  we  present  conditions  under  which  the  problem  reduces  to  the  nonnegative  inverse 
eigenvalue  problem. 


2.  The  General  Spectrum 

In  the  case  of  the  real  spectrum,  all  of  our  solutions  used  the  fact  that  the  diagonal 
(or  Jordan  canonical)  form  of  the  matrix  A  is  essentially  nonnegative,  since  all  of  its 
off-diagonal  entries  are  zero  (or  1).  When  a  matrix  C  with  real  spectrum  is  essentially 
nonnegative,  then  there  always  exists  an  r  such  that  C  +  rl  is  nonnegative.  However, 
when  A  has  complex  eigenvalues,  it  is  possible  that  it  is  never  similar  to  an  essentially 
nonnegative  matrix.  For  example,  if  the  spectrum  of  A,  u(A)  =  {1,  i,  —  ?’},  then  no  matrix 
B  exists  such  that  C  =  B~1AB  is  essentially  nonnegative.  This  is  because  there  is  no 
real  r  for  which  the  shifted  spectrum  <j(C  +  rl)  =  {1  +  r,  i  +  r,  —  *  +  r}  is  the  spectrum 
of  a  nonnegative  matrix.  (The  spectrum  a  =  {1  +  e,i,  —i}  is  in  fact  the  spectrum  of  a 
nonnegative  matrix  for  any  e  >  0.) 

It  is  clear  that  in  order  for  a  matrix  A  to  be  quasimonotone  nondecreasing  with 
respect  to  a  nonnegative  cone,  it  is  necessary  that,  for  some  real  r,  the  shifted  spectrum 
a(A  +  rl)  be  the  spectrum  of  some  nonnegative  matrix,  i.e.,  the  shifted  spectrum  of  A 
must  solve  the  nonnegative  inverse  eigenvalue  problem.  Inverse  eigenvalue  problems  are 
well-studied  problems  in  linear  algebra  (see  [5])  and  perhaps  the  most  important  unsolved 
problem  of  this  type  is  the  nonnegative  inverse  eigenvalue  problem,  which  asks  when  a 
set  of  numbers  is  the  spectrum  of  a  nonnegative  matrix.  Although  the  problem  has  been 
studied  extensively,  in  general  it  remains  unsolved  (see,  for  example,  [4]). 

We  offer  a  solution  to  our  problem  when  the  shifted  spectrum  of  A  (for  some  real 
r)  solves  the  irreducible  nonnegative  inverse  eigenvalue  problem,  in  that  cr (A  +  rl)  is 
the  spectrum  of  an  irreducible  nonnegative  matrix  An.  (The  case  where  the  shifted 
spectrum  of  A  can  be  decomposed  into  subsets,  each  of  which  solves  the  irreducible 
nonnegative  inverse  eigenvalue  problem  individually,  requires  that  we  extend  the  theory  of 
this  paper  using  techniques  similar  to  those  shown  in  [1]  for  reducible  matrices.  However, 
unfortunately  no  general  theory  exists  based  strictly  on  the  spectrum  of  A  for  such 
matrices.) 

If  there  exists  a  B  >  0  such  that  Ae  =  B~1AB  is  essentially  nonnegative  and  ir¬ 
reducible,  then  Perron-Frobenius  theory  tells  us  that  Ae  has  a  real  eigenvalue  Ai  with 
Ai  >  Re(A,)  for  i  =  2, ..., n  and  an  associated  positive  eigenvector  yi  (since  An  =  Ae+rl 
is  nonnegative  and  irreducible).  Therefore,  it  is  necessary  for  A  to  have  an  eigenvalue  Ai 
with  Ai  >  Re(Ai)  for  i  =  2,  ...,n  and  an  associated  positive  eigenvector  xi  since  xi  =  B yi 
with  yi  >  0,  B  >  0,  and  B  nonsingular.  If  A  and  Ae  are  similar,  we  show  that  it  is 
sufficient  that  A  have  a  positive  first  eigenvector  in  order  for  the  change-of-basis  matrix 
B  to  be  nonnegative,  so  that  A  is  quasimonotone  with  respect  to  a  nonnegative  cone. 
We  note  the  similarity  between  this  result  and  the  result  in  [1]  for  matrices  with  real 
spectra,  where  it  was  also  sufficient  that  A  have  a  first  eigenvector  xi  >  0. 
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3.  A  Theorem  For  The  General  Spectrum 

Theorem.  Let  A  £  Rna:n  for  n  >  2.  In  order  for  A  to  be  quasimonotone  nondecreasing 
with  respect  to  a  nonnegative  cone,  it  is  necessary  that  for  some  real  r  the  shifted  spectrum 
a(A  +  rl)  solve  the  nonnegative  inverse  eigenvalue  problem,  and  that  A  have  a  first 
eigenvalue  \\  where  Ai  >  Re(Aj)  for  i  =  2,...,n  with  an  associated  nonnegative  first 
eigenvector  xi  >  0.  Furthermore,  it  is  sufficient  that  A  be  similar  to  an  irreducible 
essentially  nonnegative  matrix  Ae  and  that.  A  have  a  positive  first  eigenvector  xi  >  0. 

Proof.  The  necessary  conditions  follow  from  the  previous  discussion,  and  we  note  that 
for  matrices  with  real  spectra,  the  requirement  on  the  shifted  spectrum  is  always  trivially 
satisfied.  To  show  the  sufficient  conditions,  let  Ae  be  similar  to  A  where  Ae  is  irreducible 
and  essentially  nonnegative,  so  that  a(A  +  rl )  solves  the  irreducible  nonnegative  inverse 
eigenvalue  problem.  Furthermore,  let  A  have  a  positive  first  eigenvector  xi  >  0  associated 
with  the  real  eigenvalue  of  greatest  real  part.  Let  D  be  the  real,  block-diagonal  canonical 
form  of  A  and  Ae  so  that  D  =  X~x  AX  =  Y~1AeY  where  the  real  and  imaginary 
components  of  the  (generalized)  eigenvectors  of  A  and  Ae  are  the  columns  of  X  and  Y 
respectively  (further,  let  the  first  eigenvector  xi  >  0  be  the  first  column  of  X).  Then 
Ae  —  YX~1AXY~1 ,  but  it  is  not  necessarily  true  that  XY~x  >  0,  as  required. 

Let  An  =  Ae  +  rl  with  r  chosen  large  enough  that  all  diagonal  elements  of  An  are 
positive  and  An  is  nonsingular.  Then  An  is  nonnegative  and  primitive  (irreducible,  with 
only  one  eigenvalue  of  maximum  modulus;  see,  e.g.,  [11]).  Then  a  further  change  of 
basis  using  (1/gi  An)k,  where  =  X\  +  r  is  the  first  eigenvalue  of  An,  yields  Ae  = 
(1/ /lzi  An)~kYX~1AXY~1(\/ gi  An)k .  Since  An  is  primitive  and  nonnegative,  then  as 
k  — >  oo,  (1/gi  An)k  — >  yizf,  where  yi  is  the  first  eigenvector  of  An  (and  of  Ae),  Zi 
is  the  corresponding  first  left  eigenvector,  yi  >  0,  zi  >  0,  and  yfzi  =  1.  Hence, 
XY-'il/^An)*  ->XlzT  >  0.  Therefore,  for  some  finite  k  this  change  of  basis  B  will 
be  nonnegative,  Ae  =  B~1AB  with  B  >  0,  and  A  is  quasimonotone  nondecreasing  with 
respect  to  a  nonnegative  cone.  □ 

We  note  that  the  above  theorem  has,  as  a  special  case,  our  previous  result  for  ma¬ 
trices  with  real  spectra,  since  such  shifted  spectra  always  solve  the  nonnegative  inverse 
eigenvalue  problem,  and  we  showed  (using  a  different  technique)  that  a  positive  first 
eigenvector  was  a  sufficient  condition  for  quasimonotonicity  with  respect  to  a  nonnega¬ 
tive  cone.  Furthermore,  our  previous  result  held  only  for  n  >  3,  as  does  this  one  as  well, 
since  a  necessary  condition  for  quasimonotonicity  (a  real  eigenvalue)  and  the  general 
spectrum  (implying  a  pair  of  complex  eigenvalues)  require  at  least  three  eigenvalues. 
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